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We present a complete derivation of two-particle states of the one-dimensional extended Bose- 
Hubbard model involving attractive or repulsive on-site and nearest-neighbour interactions. We find 
that this system possesses scattering resonances and two families of energy-dependent interaction- 
bound states which are not present in the Hubbard model with the on-site interaction alone. 
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Among the tight-binding lattice models of condensed 
matter physics, the Hubbard model plays a fundamental 
role [H- In its basic form, the Hubbard model describes 
particle tunneling between adjacent lattice sites as well 
as short range (contact) interaction between the particles 
on the same lattice site. Despite apparent simplicity, this 
model is very rich in significance and implications for the 
many body physics on a lattice . This is perhaps most 
profoundly manifested with numerous important exper- 
imental and theoretical achievements with cold neutral 
atoms trapped in optical lattice potentials 0, [1] , wherein 
the Hubbard model is being realized with unprecedented 
accuracy. 

The next level of generalization pertaining to, e.g., 
electrons in a crystal [H] or dipolar atoms 0] or molecules 
in an optical lattice, yields the extended Hubbard 
model involving longer-range interactions between the 
particles on the neighboring lattice sites. Under certain 
conditions, namely, for hard-core bosons or strongly in- 
teracting fermions at half-filling, the extended Hubbard 
model can be mapped onto various lattice spin models 
i- 

A remarkable Hubbard model effect demonstrated in 
a seminal experiment of Winkler et al. Q is the bind- 
ing of pairs of particles into composite objects by the 
on-site interaction, which can be either repulsive or at- 
tractive [loL 11 1 . In the present paper, we study the 



one-dimensional extended Bose-Hubbard model involv- 
ing attractive or repulsive on-site and nearest-neighbour 
interactions. We present a complete derivation of two- 
particle states of the system and show that it possesses 
scattering resonances and energy-dependent interaction- 
bound states which are not present in the Hubbard model 
with the on-site interaction alone. 

The Hamiltonian of the extended Hubbard model reads 



(1) 



where {bj) is the creation (annihilation) operator and 

fij = b]jbj the number operator for a boson at jth lattice 
site with energy Sj, J(> 0) is the tunnel coupling between 



adjacent sites, and U and V are, respectively, the on-site 
and nearest-neighbor interactions which can be attractive 
or repulsive. 

We seek two-particle eigenstates of Hamiltonian ([T]) 
in a homogeneous lattice, Sj = e for all j. For con- 
venience, we set e = 0, which amounts to a shift of 
the zero of energy. We then expand the state vec- 
tor in terms of the non-symmetrized coordinate basis 
{ \X]:Vr)} as I*) = Ejj' *(a;i,%') \xi,yj'), where Xj = 
dj and yji = dj' are the particle positions and d is the lat- 
tice constant. The standard (symmetrized) bosonic basis 
is recovered via the transformation \2j) = \xj,yj) and 
Ij'lj') = -^il^j^Vr) + IVj^^j')) U ^f)- Defining the 
center of mass R = ^{x + y) and relative r = x — y coor- 
dinates, the two-particle wavefunction can be factorized 
as '^{x,y) = e^^ ^ ij; K (r) , where the relative coordinate 
wavefunction t/'k (r) depends on the center-of-mass quasi- 
momentum K g [— Tr/d, Tr/d] as a continuous parameter. 
The eigenvalue problem H \'$) = E\'i>) then reduces to 
the three-term difference equation 

-JK[ipKiri-i) + ipK{ri+i)] 

+ [USr^0 + Vi6r,d + Sr,^d) ~ EK]^pK{n) =0, (2) 

with Jk = 2J cos{Kd/2) and r.; = di {i — j ~ j'). 
The above equation admits two kinds of solutions, corre- 
sponding to the scattering states of asymptotically-free 
particles and to the interaction-bound states of particle 
pairs. 

We begin with the analysis of scattering solutions 
of Eq. ([2]). They are most straightforwardly obtained 
with the standard symmetrized ansatz ipK,k{'fi 7^ 0) cx: 
g-ife|T-i| _|_ g22(5ic.fcgifc|ri| representing plane waves undergo- 
ing a scattering phase shift SK,k ■ This immediately yields 
the eigenenergies Ex^k = — 2 J/^ cos (fed), which are equal 
to the sum of Bloch bands ex,y = — 2 Jcos[(i4'/2 ± k)d] 
of two (asymptotically) free particles x, y with relative 
quasimomentum k [lol . [llj. For a given value of the 
center-of-mass momentum K , and thereby Jk, the lowest 
Exfi — and highest Ek,-k — "^Jk energy states are 

attained, respectively, at k —>■ and k — *■ n/d. The con- 
tinuum of energies EK,k and the corresponding density 
of states p{E,K) oc d(kd)/dE ^ [(2Jk)2 - £;2]-i/2 are 
shown in Fig. [T] The wavefunction of scattering states is 
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FIG. 1: Energies versus the center-of-mass momentum K for a 
pair of bosons in a ID lattice described by the extended Hub- 
bard model. The continuum spectrum corresponds to energies 
EK,k of the scattering states, with the shading proportional 
to the density of states p{E,K). The lines correspond to en- 
ergies Ek of the two particle bound states for various values 
of interaction strengths [U/ J,V/ J]. 



given by 



i'KMO) = cos{5'') — ^7^, (3a) 



cos{kd + S 



(0) ^ 
K,k) 



i^K,k{rt 7^ 0) = cos(fc|rj| + 5K,k), (3b) 
where the phase shifts 5'^^\ and 5K,k are defined through 

(4a) 



tan(4°\.) 



U 



tan((5x,fc) = 



2Jk sm{kd) ' 
JkU + [2JKC0s{kd) 



U]Vcos{kd) 



{UV - 2Jk[Jk - V cos (kd)] } sm{kd) 



(4b) 



Note that when V = 0, we have Sx.k = ^k\^ 
and the above expressions reduce to those of with 
all the consequences discussed there. Here we ex- 
amine the role of nearest-neighbour interaction V ^ 
due to which the scattering amplitude f{SK.,k) = 
if^^2iSK,k _ 2) or, for that matter, the cross-section 

o'K.k = \f {^K,k)\^ = sin^((5x.fe) can exhibit novel in- 
triguing features seen in Fig. [51 The cross-section van- 
ishes at 6K,k = 0, which, according to Eq. (|4b[) . requires 

EK.k = ic/[i ± 

cos(fco?) 



or, more explicitly, 



-U/{■iJK)[l±^/l - 8J'^/{UV)l with the con- 
dition 8J|-/(C/y) < 1. The cross section attains its max- 
imal value <7K,k = 1 at dK.k = i7r/2, which corresponds 
to a pair of impenetrable (hard-core) bosons. According 
to Eq. (|4b[) . this happens when either sin(fcd) and 
^ ^UV/{U + 2V) (for kd 0, tt, respectively), or 
EK.k = U — 2J]^/V, the last equality yielding the explicit 
condition cos{kd) — Jk/V — U/{2Jk)- 





1 -1 -0.5 

Kd/n 



FIG. 2: Two-particle scattering cross-section GK,k versus the 
center-of-mass K and relative k quasi-momenta for several 
values of interaction strengths \U / J,V / J]. Flipping simulta- 
neously the sign of both U and V is equivalent to shifting 
k ^ k + n / d. 



At the edges of the scattering band, kd 0, tt, we 
can define generalized ID scattering lengths a^^^^ 



A0,7r) 



- lini 

kd^ ' 



dSK,k ^ UV - 2Jk{Jk T V) 
dk UV ±Jk{U + 2V) 



(5) 



Thus, the scattering length uk vanishes at the bottom of 
scattering band, fc ^ 0, when Jk = \[V±yJV{2U + V)], 
and at the top of the band, k — ^ Tr/d, when Jk = 
^V{2U + V)], with the condition V{2U + V)>Q. The 
scattering length diverges when Jk = T-W for kd — > 
0,7r, where W = UV/{U + 2V). In other words, the 
divergence of a^^'^^ occurs when the center-of-mass quasi- 
momentum |iir| is equal to 



Kb. = —, arccos 
d 



W_ 
2J 



(6) 



with the condition < ^T4^/2J < 1 for kd 0,7r, re- 
spectively. As will become apparent from the proceeding 
discussion, Kn indicates the emergence of scattering res- 
onances associated with the bound states (see Fig. (H]). 

We now consider the two-particle bound states. Using 
the exponential ansatz ipKi^i 7^ 0) c< a'^ ^ in Eq. 
after little algebra, we obtain 

JKVal, + {VU - Jl)a], + Jk{V + U)aK + = 0. (7) 

Solutions of Eq. ([7]) with \aK\ < 1 determine the normal- 
ized bound states 

2JKCtK 



V'k(0)=A/' 



UaK 
^Kin^O)=Ma 

where A/" = [iPK + 2/(1 
given by 



JKia^ 
1*1-1 

K ' 



K 



1) 



-.Mlpk, (8a) 
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(8b) 

with the energies 
(9) 
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FIG. 3: The U, V diagram for the existence of the second 
bound solution |q^'| < 1 of Eq. © with Kc < \K\ < Tv/d, 
with shading proportional to Kc G [0,Tv/d]. 



Note that complex solutions of Eq. ([7|) do not correspond 
to the bound states, even if \aK\ < 1, since the energy of 
Eq. (O should be real. 

It can be shown that Eq. ([7]) admits at most two so- 
lutions corresponding to the bound states. Obviously, 
for noninteracting particles U — V = 0, there can be no 
bound state. For J7 ^ and ^ = [i3| and for [/ = 
and V 0, there is only one bound solution a^^^ at any 
K. For any other values oi U,V ^ 0, the first bound 
solution a^'' exists at any K, and the second bound so- 

(2) 

lution a)^ exists at \K\ > Kc, where the critical Kc is 
shown in Fig. [3] and is defined as 



Kc 



Kr if \W/2J\ < 1 
0^ otherwise 



(10) 



We thus see that the scattering length diverges when 
the second bound state approaches the edge of the scat- 
tering continuum. This signifies the appearance of the 
scattering resonance at |-?sr| = Kr which is determined 
by Eq. ([6|) under the condition \W/2J\ < 1. When, 
however, this condition is not satisfied, no scattering 
resonance is present, and the second bound state ex- 
ists for all K G [— vr/d, 7r/d] with the energy below or 
above the scattering continuum depending on whether 
W = UV/ {U + 2V) is negative or positive, respectively 
(see Fig.ll]). 

In general, analytic expressions for the bound solutions 
of Eq. ([7]) are too cumbersome for detailed inspection, 
but several special cases yield simple instructive results, 
(i) With only the on-site interaction U ^ Q and V = 
there is one bound solution a"^ = (C/ — Ek) / {2Jk) with 
the correspondin g en ergy Ek = sgn{U)y/U'^ + 4J^, as 
was discussed in [10|. (ii) With very strong on-site in- 
teraction |J7| ^ oo and V 0, the first bound solution 



FIG. 4: Relative coordinate wavefunctions ipKiri) of the 
two-particle bound states for several values of interaction 
strengths [[//J, V^/J]. Thick (blue) bars correspond to the 
first bound solution of Eq. ^ a.t K = Q. Thinner (red) bars 
correspond to the second bound solution at = 0, if it exists 
for all K, or a.t K — 37r/4d, if it exists for |A"| > Kc (cf. 
Figs. [U and EJ. 



finitely bound pair. The second more relevant solution 
a^' = —Jk/V and Ek = V + Jj^/V describes a pair 
of hard-core bosons, "0^(0) — 0, which are bound by the 
nearest-neighbour interaction V provided \Jk/V\ < 1 
14|. Note that in this limit we have W = V, and the 
last condition for the existence of the second bound state 
again reduces to \K\ > Kc- (iii) Finally, a rather curi- 
ous and simple case is realised with U = —V. The first 
bound solution is similar to that in (i), but with U re- 
placed by V, = (y — Ek)/{2Jk) and the energy 
Ek = sgn(y)^y^ 4- 4J|-, which corresponds to binding 
mainly by the off-site interaction. The second bound so- 



is trivial, a 



and Ek = U , representing an in- 



lution is similar to that in (ii), but now with V replaced 
by [/, a^' = -Jk/U and Ek ^ U + JfJU, provided 
\Jk/U\ < 1 (note that now W = U). This solution cor- 
responds to the on-site interaction binding. 

More generally, when the on-site U and off-site V inter- 
actions have different sign, the first bound state is asso- 
ciated mainly with the stronger (in absolute value) inter- 
action, and the second bound state with the weaker one. 
When, however, U and V have the same sign and com- 
parable strength, the bound states have mixed nature in 
the sense that both interactions significantly contribute 
to the binding. To illustrate the foregoing discussion, 
in Fig. [4] we show the wavefunctions of bound states for 
several cases pertaining to the on-site, off-site and mixed 
binding, while the corresponding energy dispersion rela- 
tions are plotted in Fig. [TJ 

To conclude, we have derived a complete solution of 
the two-body problem in a one-dimensional extended 
Hubbard model. We have found that depending on the 
strength of the on-site and nearest-neighbour interac- 
tions, this system possesses one or two families of bound 
states as well as scattering resonances corresponding to 
the degeneracy of the bound and scattering states. Our 
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results pertaining to the on- and ofF-site pairing mecha- 
nisms might be relevant to the studies of high- Tc super- 
conductivity Q and to the more recent quest to realize 
related physics with cold trapped atoms or molecules in 
optical lattices [H, [l^ . 
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